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ABSTRACT 

We study in this paper a theory of free anyons associated to free conformal field theo- 
ries defined on Riemann surfaces with a discrete and nonabelian group of authomorphisms. 
The particles are exchanged according to a nonabelian statistics, in which the R— matrix 
satisfy a multiparametric generalization of the usual Yang— Baxter equations. 
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1. INTRODUCTION 



In two recent papers, [Q and 0, we have introduced a method for constructing con- 
formal blocks with given monodromy properties, provided the monodromy properties can 
be related to those of an affine algebraic curve defined on the complex sphere CPi. The 
method was applied, as a first step, to the case in which the algebraic curves have a discrete 
group of symmetry D m . Surely applications to more complicated curves are possible, but 
we do not have the proof that the method works in the most general case. Moreover, it is 
difficult to relate our conformal blocks to a confer mal field theory at genus zero, but it is 
relatively simple to associate them with the theory of the b — c systems on algebraic curves 
with D m group of symmetry. Also in the simple case discussed, the conformal blocks, 
builded in terms of free fields and twist fields 0, that are explicitly known, exhibit 
nontrivial properties when two of the twist fields are exchanged. The twist fields turn out 
to be anions satisfying a nonabelian statistics. The purpose of this paper is to study the 
braid group statistics of the twist fields. 



2. CONFORMAL FIELD THEORIES WITH NONABELIAN GROUP OF 
SYMMETRY 



In order to fix the notation, we briefly review the conformal field theories with non- 
abelian group of symmetries introduced in and refs. M, M. We begin considering the 
case in which the nonabelian group of symmetry is simply D m . Further generalizations 
are possible as we will discuss in the next Section. 

First of all let us introduce the b — c systems of spin A (A G Z) on an algebraic curve 

S= [ d 2 zbdc + c.c. (2.1) 

"Eg is defined by the vanishing of the Weierstrass polynomial 

y 2m - 2q(z)y m + q 2 (z) - p(z) = (2.2) 

where q(z) and p(z) are polynomials of degree mr and 2r' respectively. The point z = oo 
is reached performing the substitution z' = 1/z, so that the variables z and z' represent a 
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system of local coordinates describing the complex sphere CPi. The genus g of the curve 
"Eg is given by the Riemann-Hurwitz formula: 



2g - 2 = 2m((m - l)r - 2) + 2mr' mr > r' (2.3) 

2g — 2 = 2r'(m — 1) + 2mr' — 4m mr < r' (2-4) 

The solution of eq. (|2.2[) is a multivalued function (on CPi) y^(z). The index / = 



0, . . . , 2m — 1 denotes the branches of y(z), which are defined in the following way: 




e 2 ^ y q (z) + ^(z) 0</<m-l 

2 l I /— — V / 

e~^r W n(z) — J viz] m < I < 2m — 5 



where e m = 1. The fact that the only two coordinates on E 3 are z and y(z), allows us to 
treat the b — c systems as a multivalued field theory on the complex sphere. The function 
y(z) has two different kinds of branch points. The branch points cti, i = 1, . . . , A^ a = 
max(2mr, 2r') are the projections on CPi of the zeros of y(z). The branch points flj, 
j = 1, . . . , Np = 2r' represent instead the roots of the polynomial p(z). For simplicity the 
integers r and r' are chosen in such a way that z = oo is never a branch point. The local 
group of monodromy, describing the sequence in which the branches of the curve y(z) are 
exchanged at the branch points, contains the nonabelian group D subgroup. 

If we put A > 0, a basis of meromorphic A-differentials and 1 — A-differentials in which 
all the elements have an independent behavior at the branch points, marked by the index 
< k < 2m - 1, is given by MM: 



Bf{z)dz x = (y {l) {z)y qk ' ai (p(z)) qk ^ dz x 
Ci l) (z)dz 1 - X = (y il) (z)y mqk ' ai (p(z)y qk ^ dz x ~ x 



(2.6) 



where 



and 



[fc]m + A(l 



m 



Qk, ai = [k m ] = [k + m] m = k (2.7) 

m 



A 
2 

1 - A 



qk,p s = - g k = 0, . . . , m - 1 

k = m, . . . , 2m — 1 



(21 
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In eq. ( p.6| ) all the possible monodromy properties around the branch points that a mero- 
morphic tensor can exhibit on £„ are represented. Moreover any meromorphic tensor on 
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T, g can be expanded in terms of the elements of the basis ( |2.6| ), the coefficients of this 
expansion being rational singlevalued functions of z |§. As a consequence the mpoint 
functions of the b — c systems on an algebraic curve, explicitly computed using the method 
of fermionic construction [0,0, are also decomposable in n-tensor products of the modes 
B^\z)dz x and C^\z)dz 1 ~ x . In fact, the building blocks entering the n-point functions in 
the case A > 1 are the following tensors: 

-i r \ 1 1 —\ 2m— 1 

Kf >(«,»)*»*»'-» = -L!!i*H_ £ BfMCfV) (2.9) 

/nm ? — in < * 



2m z — w 

k=0 



and the zero modes 



Q ik>k (z)dz x = z t - 1 B k (z)dz x l<i k <N bk (2.10) 
K^ 1 \z : w)dz x dw 1 ~ x is characterized by a single pole in the variable z when z = w and 



1 = 1'. In eq. ( gTTo|) the index ik runs over all the zero modes Qi kj k(z)dz x sharing the 
same monodromy properties of Bk(z)dz x . The number of these zero modes Nf, k can be 
explicitly calculated in two ways: by explicit construction of the zero modes as in ref. 
or from the residues of the ghost current j' l \z) =: b(z)c(z) : as we will see below. The 
case A = 1 can be treated in an analogous way provided we replace K^ 1 \z,w)dz x dw 1 ~ x 
with the differential of the third kind: 

V) {z)dz = Kx =^ a)dz - h)dz (2 - n) 

This differential has two simple poles in z = a and z = b when I = V and / = I" respectively. 
Moreover, when A = 1, one has to consider also the zero mode in the c— fields, given 
by Cq(z) = const. dz°, so that N Ck = 5k,o8i-\fi. This decomposition in 2m k— sectors, 
corresponding to the different elements of the basis (|2.6| ), is also evident in the vacuum 
expectation values (vev) of the ghost current and of the energy momentum tensor T(z). In 
fact, it is possible to pick up in these vev's the contributions < T^\z) > and < (z) > 
coming from the modes Bk(z) that correspond to a given monodromy behavior at the 
branch points: 

2m-l . 2m-l 



< Til) ^) >= h E < > < >=iE< J i l) w > 

(2.12) 



2m ^ h 2m 

k=0 fc=0 
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Since it will become useful for later purposes, we give the explicit expression of the vev 

<4 l \z) >: 

< J®{z) >dz = ^-d z \ogC«\z)dz (2.13) 

The vev < Tjp(z) > has been already derived in and we do not report it here. It is 
possible to consider < T^\z) > and < J k (z) > as the vev's of the energy momentum 
and of the ghost current of the fields belonging to a given k— sector, distinguished by their 
monodromy behavior at the branch points. For example the particle content of the k— th 
sector is provided by the residues of < J k \z) >. In Section 5 of ref. ]2| we showed that 
the total residue of < ■/£ (z) > at z = oo is 1 — 2A. Moreover, in the case of a curve 
with nonabelian monodromy group, there are also the residues qk, ai (I) and qk,Pi (0 at the 
branch points, which are dependent on the branch I of Jk(z) in which they are computed: 

q kiai (l)=0 0</<m-l 

(2.14) 

=Qk, ai m < I < 2m - 1 

and 

q k ,p 3 (I) = q k ,p 3 < I < 2m - 1 (2.15) 



The charges qk, ai and qk,0 t were already given in eqs. ( |2.7[ ) and ( |2.8| ). Summing the residues 



of < jjp (z) > over all the possible branch points and over the branches / = 0, . . . , 2m — 1, 
we get the following equation determining the numbers of zero modes Ni, k , N Ck in a given 
k— sector: 

2m- 1 N a Np 

Nhk -N Ck = i-2\- E^.^w-E^ ( 2 - 16 ) 

1=0 i=l j=l 

The term 1 — 2A comes from the residue of the current < jj. ! (z) > at z = oo and it is 
the usual residue one would expect in the flat case. The residues at the branch points, 
proportional to qk, ai (0 and qk,/3j , represent a nontrivial correction introduced by the global 
topology of the Riemann surface E 3 . Finally, from eq. ( |2.16| ), it is clear that the total ghost 
charge in a given k— sector requires the introduction of Nb k — N Ck zero modes, otherwise 
the amplitudes of the b — c systems vanish identically. These Nb k + N Ck zero modes are 
exactly those of eq. (|2.10| ) having the same behavior at the branch points of a given element 
B k (z)dz x of the basis (||J). 

Finally, a study of the leading order poles of < T k (z) > indicates that the effect of 
the branch points in the correlation functions can be simulated in a way that we will show 
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later by the introduction of the socalled twist fields. These fields are conformal operators 
depending on the branch points and having fractional ghost charges. In the case in which 
the monodromy group of the algebraic curve is nonabelian, the dependence of the twist 
fields on the branch points becomes necessarily nonlocal. 

One possibility to construct effectively ab — c theory on S 3 in which the splitting of the 
fields is explicit, consists in expanding the fields in powers of z and y(z). This is the most 
general expansion allowed on an algebraic curve. Moreover, on an algebraic curve ( |2.2| ) 
there are only 2m rationally independent functions fk(z), which are solution of a Riemann 
monodromy problem. The function y(z) and its powers should be linear combinations of 
the fk(z)'s, the coefficients being rational functions of z as explained above. Expanding 
also the rational functions in powers of z and collecting all the terms in fk(z) together, 
< k < 2m — 1, we get the final form of the fields: 

2m— 1 oo 

E E *"'"*/*(*)&*,< ( 2 - 17 ) 

k=0 i= — oo 
2m — 1 oo 

E E z ~ i+X ~" -fk{z)ck,i (2-18) 

A;=0 £=— oo 

Here b^,i and are arbitrary coefficients that in the quantum case become creation and 
annihilation operators of the fields propagating with the same monodromy properties of 
fk(z)- This is the scheme used in || in order to construct an operator formalism for the 
b — c systems on Riemann surfaces. Still there is an infinite number of possible solutions 
fk(z) of the Riemann monodromy problem, which are equivalent up to rational functions 
of z. The proper solutions for the b — c systems are given in ( ^.6|) . Any other basis would 
produce an operator formalism in which the residues ( p.!4| ) and (|2.15| ) of the currents 
Jk(z) are shifted by integers. However this would be in contradiction with eqs. ( |2.14|) , 
( |2.15|) obtained computing explicitly the vev < Jk(z) > with the method of the fermionic 
construction 0. 

3. BRAID GROUP STATISTICS OF THE b-c SYSTEMS 

Another possibility of realizing the theory of the b — c systems on an algebraic curve 
in agreement with eq. ( |2.12| ) and ( |2.16| ) is based on bosonization The details are 

described in ref. . See also , |T(| for a study of the Z n symmetric curves. Here we will 
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mainly study the braid group statistics of the twist fields. First of all we associate to each 
fc-sector with different monodromy properties at the branch points the couple of free b — c 
fields bk(z)dz x and c k (z)dz 1 ~ x . These fields take their values on the complex sphere and 
act on the vacua |0 > k: < k < 2m — 1. The fields belonging to different k— sectors do not 
interact. The presence of the branch points in the amplitudes is simulated by the twist 
fields (a.i) and V k (fy) for each sector k. These are primary fields as it is shown in 
ll and 0. The two point function of the fields bk and c k is defined as follows: 



Gf k \z,w)dz x dw 



1-A 



B k\ Z ) C k \ W ) TT ( Z ~ Z k ,s) TTf (W - W k ,s) , A 



Nr. 



Z — W 



n 



n 



L (w - z k , s ) (z - w kyS ) 



-dz dw 



1-A 



(3.1) 



The meaning of eq. (|3.1| ) as a two point function is due to the fact that this meromorphic 
tensor yields the vev's of the currents J k (z) and of the energy momentum tensors T k (z) 
discussed in the previous Section (see also 0). Strictly speaking, however, this is not a 
true Green function on the algebraic curve because the structure of the poles is wrong. 
For example the pole in z = w does not occur only when / = I' as it should be. However 
it is a Green function on the complex plane for the multivalued sector k. 

It was shown in that the explicit expression of eq. ( |3.1|) in terms of twist fields and 
free b — c fields is given by: 



(z, w)dz x dw 



N bk N Ck Na Np 

< o\b k (z)c k (w) n Kz kia ) n c{z k>a >) n vm n ^(^)io > k 



1-A 



s = l 



s' = l 



1=1 



3 = 1 



k<o\u b( Zk , s ) n <z k ,s') n v k {oi) n nc^oio > k 

s=l s'=l i=l j=l 



(3.2) 



Bosonizing the fields b k (z) and c k (z) with the rules: 



b k (z) = e-'^M c k (z) = e lVk{z) 

we get the following form of the twist fields: 

,({) 



< ip k (z)(p kl (z') >= -5 kk >log(z - z') 

(3.3) 



K ■ («») = ex P 



dt < J ( k l \t) > <p k {t) 



Co 



(3.4) 



V k il \pj) = exp 



i.<D dt < 4° (t) ><p k {t) 

Op* 



(3.5) 
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The vev of J^\ z ) is given explicitly in eq. ( |2.13j ) and C Qi , Cp. represent small paths 
surrounding the branch points CKj and /3j respectively. The operators appearing in eqs. 
and ( p.5| ) satisfy a nonstandard statistics ||: 

vt l \i)vt l '\i) = 

1 vPMvPi-y) 



exp 



-<?/c,7<?/c, 7 ' f ds < Jf?\s) > f ds' < jj! \s') > log(s - s') 

(3.6) 

where 7 and 7' take their values on the set of branch points {di,j3j}. In eq. (|3.6|) , the 
dependence in the statistics of the twist fields on the branch points is reminiscent of the non- 
flat spacetime background on which the theory was originally defined. The (/-parameter is a 
complicated convolution of three distributions. This is not in contrast with ref. || , where it 
was shown that the dependence of the q— parameter on the spatial coordinates is excluded 
in the flat case due to the translational invariance of the local exchange relations. First 
of all we are here on a Riemann surface with a nontrivial geometry. Secondly, the algebra 
(PUD, defined on a local patch of the Riemann surface which is topologically equivalent to 
a complex plane with punctures at the branch points 7 and 7', becomes: 

^ ( °(7)vf V) = j^W^'WvpMvPi-y) (3.7) 



with 7 7^ 7' and I ^ V . To derive eq. (|3.7|) we have simply evaluated the two line integrals 
occurring in eq. ( |3.6| ) using the following two equations, that express the coefficients 
appearing in eq. (|2.2|) in terms of the branch points: 



q\z) - p(z) = l[(z - a,) p(z) = ]J(z - ft) (3.* 

i=l i=j 




\*=i 3=1 

The remaining task is just a computation of residues. As a result, we see that the local 
algebra ( |3.7| ) has a q— parameter exp [iiTqk^ityqk,^' (I')] in which the dependence on the 
spatial coordinates is only appearing in the indices of the charges <?fc, 7 (0 and g^y^'). 



At this point it is convenient to forget for the moment the index k of the monodromy 
sector and to denote all the branch points with the symbol 7$, 1 < i < N a + Np, in such 
a way that: 

7i =OLi 1 < i < N a 

(3.10) 

=Pi-N a N a + l<i<N a + Np 

Moreover we group the the two indices U and 7$, describing the branch of V fe (7^) and 
the branch point 7$ simulated by the twist field respectively, in an unique index. The new 
composed indices will be denoted using the capital latin letters I, J,K . . .. In this condensed 
notation we have for example that (7*) = V 1 , (77) = V J , 1 < i,j < N a + Np. 
Analogously qk ni (h) = qi and so on. A sum over the index i", containing the indices U and 
7»j 1 < i < A'q, + should be understood as a sum over the two indices i and l{. These 
two indices are totally independent. The subscript i in Zj is just a convenient notation 
when many twist fields are considered altogether, otherwise we could set U = a, lj = b and 
so on, a, b = 0, ... , 2m — 1. The compact form of eq. ( p.7|) becomes: 



ylyJ = e inqiqjyJyI j^j ( 311 ) 



In eq. Q3.11|) we have still not taken into account the radial ordering (or time ordering) 



that it is understood in eq. (|3.2| ). As usual, two branch points are said to be radial ordered, 
i.e. 7 -< 7', iff I7I < |7'|, the symbol | | denoting the modulus of a complex number. Using 
this definition of time-ordering and eq. ( |3.11|) we get: 

T (V T V J ) = [e-™ qiqj 8(I -J)+ e +inqiqj 9(J - I)] T (V^ 1 ) (3.12) 

where 9(1 — J) = 9(\^i\ — is the Heaviside theta function. The algebra ( |3.12| ) is 
trivially associative despite of the fact that the (/-parameter exp(i7rqjqj) depends also on 
the indices / and J. In fact one can easily see that the following equation holds 

T (V'V'V^ 4 ) = exp [z7T (e(J, J)q iq ,j + e(J, M) qi q M + e( J, M) qj q M )] T (V M V J V T ) 

(3.13) 

independently of the way in which the twist fields in the lhs are permuted in order to get 
the rhs. In ( |3.13| ) we have defined: 

e(I, J) = + 1 if 7, -< nfj 

(3.14) 

= - 1 if 7j -< 7^ 
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The relations (|3.12| ) resemble the equations defining the quantum hyperplane of [11|, based 



on the quantum group GL q (N), where q is a root of unity. We can interpret eqs. (|3.2| ) 
and (|3.12[) in the following way. In the correlators of the b — c systems on a Riemann 



surface with nonabelian monodromy group the dependence on the branch points becomes 
no longer local, but it is expressed through the twist fields (|3.4|) and (|3.5|). These operators 
replace the coordinates on and f3j and can be considered as quantum coordinates. One is 
tempted therefore to treat the twist fields formally as coordinates of a complex hyperplane 



and to rewrite eq. (3.12) in a matricial form following refs. [1 I 



T (V T V J ) = Qmn(<1i, Qj)T {V m V n ) (3.15) 

M,N 

where 

Qmn(Qi, Qj) = exp(m?/?j)^ n [exp(-iirq kjli (h)q k)1] (lj))] 8 l fJ l { m (3.16) 

Moreover R(q), with q independent of the row indices satisfies the quantum Yang- 
Baxter equation: 

R12R23R12 = R23R12R23 (3-17) 



and the sum over M and N in eq. (|3.15| ) is intended as a sum over the small indices 



m, n 



and l m , l n . The -R mn component of Q T A j N acts on the indices i,j, m, n denoting the branch 
points 7»,7j, 7m and 7 n , while the permutation matrix 6^ acts on the branch indices 
li,lj,l m and l n . Explicitly eq. ( |3.15| ) reads: 

N a +Nfs 

m,n=l 
2m-l 

Rg n [eM-i^MQk, 7 A l j))] E S t S t T { v{lm) (^)V^( ln )) (3.18) 

where we remember that the indices l m and l m are actually independent of m and n. We 
can also generalize eq. Q3.18 ) allowing for q— parameters which are not roots of unity. In 



this case exp(i%qiqj) can be substituted by any function f(qi,qj) such that f(qi,qj) = 
f(Qj,Qi) 

The only difference between eq. ( |3.15|) and the quantum hyperplane is that the 
q— parameter appearing in the matrix Qmn depends now also on the row indices / and J. 
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The associativity property ( |3.13| ) requires that the matrix Q of eq. ( j3.16|) still satisfies a 
relation similar to the Yang-Baxter equations. Indeed it is easy to prove that: 



QMqi, q2)Q23(q2, (fe)Ql2(gi,g2) = <?23(92, 93)Ql2(?l, 92)Q23(?2, 03) (3-19) 

This is not the usual Yang-Baxter equation, but a generalization of it that was firstly found 
in ref. [12]. In our case the form of eq. ( |3.19| ) is dictated by the fact that the twist fields 



obey the nonabelian statistics ( |3.11| ), which is a consequence of the nonabelian monodromy 
group of the algebraic curve ( p.2|) . It is remarkable the fact that in refs. |13| , fl4]] . Explicit 



solutions of eq. ( |3.19| ) were found from integrable models in which the spectral parameter 
is defined on algebraic curves of genus g > 1. Multiparameter deformations of the quantum 
hyperplane such as in eq. ( |3.12| ) were also discussed in |0| and |jIB| . 



We conclude remembering that algebraic curves have many applications to knot theory 
jlTj j . Here we have expressed in terms of twist fields the braiding properties of the curve at 
the branch points. The situation in which the curve becomes degenerate, for example when 
two or more branch points coincide, is the most interesting for knot theory. However this 
situation is difficult to realize using the twist fields because the OPE's between two twist 
fields need some regularization when both of them are concentrated at the same branch 
point. 
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